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Abstract 

Conditions are studied under which there can exist a quasiequilibrium mixture of 
itinerant and localized bosonic atoms in an optical lattice, even at zero temperature 
and at integer filling factor, when such a coexistence is impossible for an equilibrium 
lattice. The consideration is based on a model having the structure of a two-band, or 
two-component, boson Hubbard Hamiltonian. The minimal value for the ratio of on- 
site repulsion to tunneling parameter, necessary for the occurrence of such a mixture, 
is found. 
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1 Introduction 



Optical lattices provide exceptional opportunity for creating various states of periodic matter 
[1-5]. In this paper, we consider N bosonic atoms in a lattice of N L sites, with the filling 
factor v = N/N L . The system is characterized by a boson Hubbard Hamiltonian, with 
on-site repulsion U, tunneling parameter J, and the number of nearest neighbors Zq. As 
is well known, Bose atoms at zero temperature and integer filling factor, form either Mott 
insulating state or superfluid state, depending on the ratio u = U/z J of the on-site repulsion 
U to the product of the tunneling parameter J and the nearest-neighbor number z . For 
a cubic three-dimensional lattice, with z = 6 and the unity filling factor [y = 1), the 
second-order phase transition between superfluid and insulating states occurs at u c = 4.9, as 
follows from strong-coupling perturbation theory [6,7] and Monte Carlo simulations [8-10]. 
At finite temperature and/or noninteger filling factor, there is the coexistence of localized 
and delocalized atoms [11-14]. 

Our concern in the present paper is to analyze a possible coexistence of delocalized, 
wandering, atoms and localized atoms in a lattice with an integer filling factor at zero tem- 
perature. As is known from the previous results, such a case cannot occur in an equilibrium 
lattice. Therefore, we need to keep in mind a kind of a nonequilibrium system. 

Let us assume that a nonequilibrium state has been prepared, where a portion of atoms 
is localized and another portion is not. This could be achieved in the process of loading 
atoms into the lattice. A nonequilibrium loading of atoms into a double-well optical lattice 
has been studied in Refs. [15-17]. Suppose that the process of such a nonequilibrium loading 
lasts the time t non that is longer than the local-equilibrium time ti oc , but shorter than the 
relaxation time t re i that is necessary for the system for passing to the total equilibrium, 

tloc "C <C t re l ■ 

In that case, in the interval of time t non <C t <C t re i, the system can be treated as quasiequi- 
librium, so that the components of the itinerant and localized atoms are in equilibrium with 
each other, while the system as a whole has not yet been equilibrated, but changes slowly. 

Note that here we consider the case of atoms inside a prescribed optical lattice. That is, 
the considered system is an artificial periodic structure, usually, of mesoscopic or nanoscopic 
size [18]. Such a setup is different from the case of a self-organized crystalline lattice of a 
quantum crystal [19-21], in which there can occur jumps of particles, connected with self- 
diffusion [22,23]. 

We leave aside the problem of how the desired quasiequilibrium structure could be cre- 
ated. Fortunately, optical lattices are highly regulated objects, whose parameters can be 
varied in a wide range [1-5]. We assume that such a quasiequilibrium system can be formed. 
But, since the quasiequilibrium has been assumed, this imposes restrictions on the system 
parameters at which the possible coexistence of itinerant and localized atoms could be real- 
ized. Our aim is to find out what are these restrictions and, in particular, what should be 
the interaction parameter u = U/Jz , when the coexistence would be admissible. 



2 Two-Band Model 

To describe the desired coexistence of atoms, we keep in mind a kind of a two-band, or 
two-component, Hubbard Hamiltonian [5], in which one band corresponds to delocalized, 
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conducting, atoms, while another band, to localized, bound, atoms. The system, as a whole, 
contains N Q + Ni delocalized atoms, where N is the number of condensed atoms and Ni, 
the number of uncondensed atoms. The number of localized atoms is N 2 . So that the total 
number of atoms is 

N = N + N 1 + N 2 . (1) 

The filling factor is the ratio 

of the total number of atoms to the number of lattice sites. The corresponding atomic 
fractions 

iV N, N 2 

satisfy the normalization condition 

n + ni + n 2 — 1 , (4) 

following from Eq. (1). 

The field operator of itinerant atoms, in order to correctly describe the Bose-condensed 
system, is represented by the Bogolubov shifted form [24] 

Cj = r) + Cj , (5) 

with the index j = 1,2, .. . , Nl enumerating lattice sites. Here, r\ is the condensate order 
parameter, defining the condensate density \i]\ 2 , and Cj is an operator of uncondensed atoms. 
Statistical averages for the operators of uncondensed atoms, Cj, and for those of localized 
atoms, bj, are such that 

( Cj ) = fa) = . (6) 
Thus, the number of itinerant condensed atoms is 

No = Yl M 2 = ™oWl • (7) 
j 

The number of itinerant uncondensed atoms is 

N, = (iVx) = u ni N L , (8) 

with the number operator 

N l = J2ch- (9) 

j 

And the number of localized atoms is 

N 2 = (N 2 ) = vn 2 N L , (10) 

with the number operator 



N 2 = J2%. (11) 
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The energy Hamiltonian has the form of a two-band Hubbard model 

U 

C l C j + 

<*J> 



h = - J E £ fe + I E ( fi 5 6 &^ + 26 5 6 ^ + ¥>a) » (12) 



where in the first term the summation is over the nearest neighbors. For simplicity, the equal 
on-site interactions are taken for all atoms. Constructing the grand Hamiltonian, we have 
to take into account the given normalization conditions, uniquely defining a representative 
ensemble for the system with broken gauge symmetry [25-27]. Then the grand Hamiltonian 
reads as 

H = H- l i N Q -n 1 N 1 -n 2 N 2 , (13) 

in which the Lagrange multipliers /io, fJ>i, and f/,2, play the role of partial chemical potentials 
guaranteeing the validity of normalizations (7), (8), and (10). The system chemical potential 
is 

fx = fi n + /iim + /i 2 n 2 . (14) 

However, in the considered case, not all these multipliers are independent. The restriction 
comes from the assumption that the system is in quasiequilibrium, such that the variation 
of a thermodynamic potential be zero. From this condition, we have [5] the relation 

M^o + n 1 ) = fi n + n x n x , 

which gives us the chemical potential 

/X ™0 + ^1™1 

\i = ■ = ii 2 ■ (15) 

n + m 

Atoms from different bands are assumed to be weakly correlated, such that 

Then the grand Hamiltonian (13) takes the form 

H = H del + H loc - vn x n 2 NU , (17) 

where the first term describes delocalized itinerant atoms, the second term, localized atoms, 
while the third term is responsible for the interactions of atoms from different bands. 

The Hamiltonian of delocalized atoms, substituting there the Bogolubov shift (5), be- 
comes the sum 

4 

H del = j2 H(n) > ( 18 ) 

n=0 

with the zero-order term 

H<® = (-zqJ + ^vn - fj^j n N , (19) 

first-order term 

H (1) = , (20) 
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second-order term 

CjCj + (2vn U + vn 2 U - Hi) c [ f j + 

(ij) 

third-order term 



H {2) = - J c] c i + ( 2 " n o u + vn 2 U - /xi) Y c ) c i + un ojYl ( C H + c i c i) > ( 21 ) 



= y^T U fofo + c h c i) > ( 22 ) 

3 

and the fourth-order term 

2 



^ (4) = ^E c kV,- (23) 



The Hamiltonian of localized atoms is 

H loc = Yl H 3' (24) 

3 

where 

H 3 = \ *»J^- (&& " l) + Hno + rn)tf - h)h 3 . (25) 

3 Itinerant Atoms 

For the operators of delocalized itinerant atoms, we can invoke the Fourier transformation 

9 = 4rE a ^ k ' aj ' (26) 



in which the summation is over the Brillouin zone and aj is a lattice vector. For concreteness, 
we consider in what follows a cubic lattice with the lattice distance a, for which 

j V v ( N 



with d being the spatial dimensionality. 

The second-order term (21) transforms to 



J7» = e 

k 

The third-order term (22) is 

U 



2 J ^ cos(k a a) + 2un U + z/n 2 t/ — /ii 



4a fc . (27) 



# (3) = jj- V^o Yl ( a l a P a k+ P + 4+ P a P a fc) • ( 28 ) 



kp 



And the fourth-order term (23) becomes 



# (4) = ^E^ww (29) 



The following consideration is in line with the general self-consistent approach to Bose 
systems with broken gauge symmetry, advanced in Refs. [25-30], and employed in Refs. 
[31-35]. We use the definition 

u k = —2 J cos(k a a) + + n + ni)U — fii (30) 

a 

and 

A = v(n + a)U , (31) 

where 

Wl = ]vE n * (n k = {a{a k )), 
k 



o = jf^2°k (<r fc = (a_ fe a fe )) . (32) 
Employing the Hartree-Fock-Bogolubov approximation yields 

Add = ^jjfb + y^u; fc 4 Q fc + \ ^ ( Aa I Q -fc + A*a_ fc a fc ) , (33) 

k k 

where the nonoperator term is 

E HFB = if (°> - iVz/ (2n? + \a\ 2 ) ^ . (34) 
Diagonalizing this Hamiltonian with the help of the Bogolubov canonical transformation 

a k = u k /3 k + v*_ k ft_ k , 

gives the Bogolubov Hamiltonian 

H del = E B + Y,£kPlPk, (35) 
k 

with the nonoperator term 

E B = E HFB + - ^( £fc - uj k ) (36) 
k 

and the excitation spectrum 

e k = sJujI - . (37) 
In equilibrium, the equation of motion for condensate atoms [5] reduces to 

dH 



w o/ »■ < 38 > 

which gives the condensate chemical potential 

pL = -z J + vXJ (no + 2ni + a) . (39) 
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The condition of condensate existence [5] requires that the spectrum be gapless, so that 



lim e k = (e k > 0) , 

which is equivalent to the Hugenholtz-Pines relation [36] and yields 

/ii = —z J + vU{l + n\ — a) . 

Introducing the notation 

e k = cos{k a a)} = 4J^sin 2 {^f\ 

transforms Eq. (30) into 

uj k = e k + A . 

The excitation spectrum (37) takes the Bogolubov form 

e k = a/ e k (e k + 2 A) . 
The normal and anomalous averages (32) are given by the integrals 

i r dk l r dk 

in which the integration is over the Brillouin zone and 

— — mth ( — —\ — — (jfc = — pot In ( — — 



n k = coth , 

2e k \2T 



coth , 
2e k \2T 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



4 Localized Atoms 

In the Hamiltonian of localized atoms (24), the site terms (25) can be represented as 



with the density operator 



U [vno + vrii- - ) - ji 



N j = b]bj . 



The average number of localized atoms per site is 



Tre-P H i 

Since the eigenvalues of Hamiltonian (46) are 



E, 



U 



— n 2 + 



U ( vn Q + vn x - - ) - /i 



n , 



(46) 
(47) 

(48) 
(49) 
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where n — 1, 2, . . ., the average number of localized atoms (48) is 



vn 2 



V°° e 



At low temperature, such that T <^ U, the maximal contribution into the 
is given by the term with n = n e ff, for which E n is minimal, that is, when 

dE n d 2 E n 
— - = > 

on on z 

The latter gives 

n eff = 2 ~ u ( n o + n i) + jj i 
under the condition U > 0. Thus, we get 

z/n 2 ~ n eff (T < [/) . 

Combining Eqs. (52) and (53) yields 

H=( v -'t\u (T < U) . 



According to Eq. (31), we have 

A 

a + n = — . 

uU 

Using this in the chemical potentials (39) and (41) results in 

/i = — zqJ + A + 2uriiU , 

Hi = -ZqJ - A + z/(l + n + ni)U 
Equating expressions (15) and (54) gives 

v — - ) U — 



2 J n + rii 

^From Eqs. (56) and (57), we find 

n (1 + 2vn 1 )U - 2A - 2z J 



n x (2u - Qvri! - 1)U -2A + 2z J ' 
The latter, employing the relation 

n 2 — 1 - n - ni , 

results in the equation 

U - 2z J + 2vU(4n 1 - 2n\ - 1) + 2A(1 - 2n x ) 
712 ~ U - 2z J + 2uU(3rn - 1) + 2A 

defining the fraction of localized atoms n 2 . 
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5 Dimensionless Equations 

It is convenient to pass to dimensionless quantities by introducing the dimensionless inter- 
action parameter 

u =h (61 > 

and the dimensionless sound velocity squared 

^h- (62) 

Also, considering a cubic lattice, we use the dimensionless wave vector with the components 

q a = -k a (a = 1,2, ...,d) . (63) 

71 

Measuring all energy quantities in units of z J, we have, instead of Eq. (30), 

u k = e k + s 2 , (64) 

where 

e k = sin2 (| 9a) » ( 65 ) 

and the Bogolubov spectrum (43) becomes 



£k = \/e k {e k + 2s 2 ) . (66) 
With the relation ga d = u, Eqs. (44) are reduced to 



X ,1 ,1 1 rl rl 



m = - ... I n k dqi . . . dq d , a = - / . . . / o k dq x ...dq d , (67) 
v Jo Jo v Jo Jo 

where and Ofc are given by Eq. (45). In particular, at zero temperature, we have 

n t = i(|-l), = <T = 0). (68) 

Equation (31) for the sound velocity becomes 

s 2 = (n + a)vu . (69) 

And for Eq. (60), we find 

2 + (2u - l)u - 4imi(2 - ni)u + 2s 2 (2nx - 1) 
U2 ~ 2 + (2v - Tju - 6un lU - 2s 2 ' ^ ' 
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6 Numerical Results 



We accomplish numerical calculations for a three-dimensional cubic lattice (d = 3, z = 6), 
with unity filing factor [y = 1), at zero temperature (T = 0). Explicitly, we solve the system 
of equations 

n i = o / ( — - 1 ) dq 1 dq 2 dq 3 , 

s 2 f 1 1 
a = - — — dqxdq 2 dq 3 , 
* Jo £ k 

uj k = e k + s 2 , e k = \/e k (e k + 2s 2 ) , 

o 3 

et= 3L sm (2 9a J ' k a = - q, - 

s 2 = (no + a)u , n — 1 - n x — n 2 , 
2 + u - Aunx(2 - m) + 2s 2 (2n 1 - 1) 
2 + u — Qun x — 2s 2 

The results for the fractions of condensed atoms, no, uncondensed atoms, n±, and localized 
atoms, ri2, are shown in Fig. 1. 

When the interaction parameter (61) is smaller than 2.75, there are no positive solutions 
for n 2 , so that the sole admissible solution is n 2 = 0, which corresponds to the single-band 
Hubbard model. The mixture can exist only for u > 2.75. This gives the lower boundary for 
the possible quasiequilibrium coexistence of itinerant Bose-condensed atoms and localized 
atoms. 



7 Conclusion 

Conditions, under which a quasiequilibrium system of coexisting itinerant and localized Bose 
atoms could be created, are analyzed. The consideration is based on a kind of a two-band, or 
two-component, boson Hubbard model [5]. One component corresponds to itinerant atoms 
and is described by the self-consistent Hartree-Fock Bogolubov approximation [25-30]. This 
approximation is known [28-30] to be well suited for superfluids. But, since it explicitly takes 
into account the global gauge symmetry breaking, it is not suitable for the Mott insulating 
state, where the gauge symmetry breaking is absent [37]. Therefore, the localized atoms 
in the insulating state are characterized, in our model, as bound atoms without tunneling 
between lattice sites. 

For numerical calculations, we take a three-dimensional cubic lattice. Lower- dimensional 
systems may have no Bose-Einstein condensate and require separate consideration [38]. 

We find that the mixture of itinerant and localized atoms can be formed only when the 
on-site interaction is sufficiently strong, such that 

u > 2.75 . 

Recall that here a quasiequilibrium system has been considered. An equilibrium system 
at zero temperature and the unity filling factor, as we know, is superfluid below u c = 4.9 
and is insulating above the latter value, undergoing, at u c , a second-order phase transition. 
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In the model, we have considered, the itinerant and localized atoms are not spatially 
separated. This distinguishes our case from that studied in Ref. [39], where superfluid 
droplets were separated in space, being surrounded by normal, nonsuperfluid, phase. 

We do not consider here the ways of creating such a quasiequilibrium mixture. Clearly, 
this should involve nonequilibrium preparation of this state. It is possible to invoke external 
temporal modulation of the system parameters. Different variants of such a modulation of 
the parameters of trapped atoms are now discussed in literature [40-45]. 

In order to better understand the physics of the mixture, let us give the following picture. 
Imagine an optical lattice characterized by the tunneling frequency J and on-site repulsion 
U . These energy parameters are related to the corresponding typical times. The energy 
ZqJ defines the time t osc = l/z^J, during which an atom oscillates in a potential well of a 
lattice site. The smaller zqJ, the longer this oscillation time. In the limit of no tunneling, 
this time is infinite. The energy U defines the wandering time t wan = 1/U, during which 
an atom realizes a jump between the neighboring lattice sites. The larger U, the shorter 
this wandering time. In the limit of infinite U, there are no jumps, all atoms are completely 
localized, and the wandering time is zero. The ratio of the oscillation time to wandering 
time is exactly the parameter u. The condition that u > 2.75 implies that 



twan 

This means that the oscillation time has to be sufficiently long as compared to the wandering 
time. In the other case, atoms could not be localized. 

Acknowledgement 

One of the authors (V.I.Y.) is grateful for financial support to the Russian Founda- 
tion for Basic Research and another author (A.R.) acknowledges financial support from the 
Volkswagen Foundation. 



11 



References 

[1] D. Jaksch and P. Zoller, Ann. Phys. (N.Y.) 315, 52 (2005). 

[2] O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 (2006). 

[3] C. Moseley, O. Fialko, and K. Ziegler, Ann. Physik 17, 561 (2008). 

[4] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008). 

[5] V.I. Yukalov, Laser Phys. 19, 1 (2009). 

[6] J.K. Freericks and H. Monien, Phys. Rev. B 53, 2691 (1996). 

[7] N. Elstner and H. Monien, Phys. Rev. B 59, 12184 (1999). 

[8] S. Wessel, F. Alet, M. Troyer, and G.G. Batrouni, Phys. Rev. A 70, 053615 (2004). 

[9] B. Capogrosso-Sansone, N. V. Prokofiev, and B. V. Svistunov, Phys. Rev. B 75, 134302 
(2007). 

[10] B. Capogrosso-Sansone, S.G. Soyler, N. Prokofiev, and B. Svistunov, Phys. Rev. A 77, 
015602 (2008). 

[11] D.B. Dickerscheid, D. van Oosten, P.J. Denteneer, and H.T.C. Stoof, Phys. Rev. A 68, 
043623 (2003). 

[12] A. E. Meyerovich, Phys. Rev. A 68, 051602 (2003). 

[13] L. I. Plimak, M. K. Olsen, and M. Fleischhauer, Phys. Rev. A 70, 013611 (2004). 

[14] F. Gerbier, Phys. Rev. Lett. 99, 120405 (2007). 

[15] VI. Yukalov and E.P Yukalova, Phys. Rev. A 78, 063610 (2008). 

[16] V.I. Yukalov and E.P. Yukalova, Laser Phys. Lett. 6, 235 (2009). 

[17] V.I. Yukalov and E.P. Yukalova, Phys. Lett. A 373, 1301 (2009). 

[18] F.A. Buot, Phys. Rep. 234, 73 (1993). 

[19] V.I. Yukalov, Ann. Physik 36, 31 (1979). 

[20] V.I. Yukalov, Ann. Physik 38, 419 (1981). 

[21] V.I. Yukalov and V.I. Zubov, Fortschr. Phys. 31, 627 (1983). 

[22] V.I. Yukalov, Physica A 89, 363 (1977). 

[23] D.I. Pushkarov, Quasiparticle Theory of Defects in Solids (World Scientific, Singapore, 
1991). 

[24] N.N. Bogolubov, Lectures on Quantum Statistics (Gordon and Breach, New York, 1970), 
Vol. 2. 

12 



[25] V.I. Yukalov, Phys. Rev. E 72, 066119 (2005). 

[26] V.I. Yukalov, Laser Phys. 16, 511 (2006). 

[27] V.I. Yukalov, Int. J. Mod. Phys. B 21, 69 (2007). 

[28] VI. Yukalov, Phys. Lett. A 359, 712 (2006). 

[29] VI. Yukalov, Laser Phys. Lett. 3, 406 (2006). 

[30] VI. Yukalov, Ann. Phys. (NY.) 323, 461 (2008). 

[31] VI. Yukalov and H. Kleinert, Phys. Rev. A 73, 063612 (2006). 

[32] V.I. Yukalov and E.P Yukalova, Phys. Rev. A 74, 063623 (2006). 

[33] V.I. Yukalov and E.P. Yukalova, Phys. Rev. A 76, 013602 (2007). 

[34] A. Rakhimov, C.K. Kim, S.H. Kim, and J.H. Yee, Phys. Rev. A 77, 033626 (2008). 

[35] A. Rakhimov, E.Y. Sherman, and C.K. Kim, Phys. Rev. B 81, 020407 (2010). 

[36] N.M. Hugenholtz and D. Pines, Phys. Rev. 116, 489 (1959). 

[37] I. Danshita and P. Naidon, Phys. Rev. A 79, 043601 (2009). 

[38] VI. Yukalov and M.D. Girardeau, Laser Phys. Lett. 2, 375 (2005). 

[39] F. Cinty, P. Jain, M. Boninsegni, A. Micheli, P. Zoller, and G. Pupillo. larXiv: 1005.24031 
(2010). 

[40] K.P Marzlin and VI. Yukalov, Eur. Phys. J. D 33, 253 (2005). 

[41] V.O. Nesterenko, A.N. Novikov, F.F. de Souza Cruz, and E.L. Lapolli, Laser Phys. 19, 
616 (2009). 

[42] LY. Yang, L.B. Fu, and J. Liu, Laser Phys. 19, 678 (2009). 

[43] V.I. Yukalov, E.P. Yukalova, and V.S. Bagnato, Laser Phys. 19, 686 (2009). 

[44] A. Zenesini, H. Lignier, C. Sias, O. Morsch, D. Ciampini, and E. Arimondo, Laser Phys. 
20, 1182 (2010). 

[45] V.I. Yukalov, Laser Phys. Lett. 7, 467 (2010). 



13 



Figure Caption 



Fig. 1. Atomic fractions, as functions of the dimensionless interaction parameter u, 
for the condensed atoms (solid line), uncondensed atoms (dotted line), and localized atoms 
(dashed-dotted line). 
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Figure 1: Atomic fractions, as functions of the dimensionless interaction parameter u, for 
the condensed atoms (solid line), uncondensed atoms (dotted line), and localized atoms 
(dashed-dotted line). 
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